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Abstract 



We investigate the large time behavior of an axisymmetric model for the 3D 
Euler equations. In [22J, Hou and Lei proposed a 3D model for the axisymmetric 
Oh ■ incompressible Euler and Navier-Stokes equations with swirl. This model shares 

many properties of the 3D incompressible Euler and Navier-Stokes equations. 
The main difference between the 3D model of Hou and Lei and the reformulated 
3D Euler and Navier-Stokes equations is that the convection term is neglected 
in the 3D model. In |24j . the authors proved that the 3D inviscid model can 
develop a finite time singularity starting from smooth initial data on a rectangular 
domain. A global well-posedness result was also proved for a class of smooth 
initial data under some smallness condition. The analysis in [24j does not apply to 
the case when the domain is axisymmetric and unbounded in the radial direction. 
In this paper, we prove that the 3D inviscid model with an appropriate Neumann- 
Robin boundary condition will develop a finite time singularity starting from 
smooth initial data in an axisymmetric domain. Moreover, we prove that the 3D 
inviscid model has globally smooth solutions for a class of large smooth initial 
data with some appropriate boundary condition. 



X ' 1 Introduction 



Whether the 3D incompressible Navier-Stokes equations can develop a finite time sin- 
gularity from smooth initial data with finite energy is one of the most challenging 
questions in nonlinear partial differential equations [13] . There have been many pre- 
vious studies devoted to this challenging question, see e.g. [U El El [HI [121 E] and two 
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recent review articles [U [5] . There have been a number of attempts to investigate pos- 
sible finite time singularity formation of the 3D Euler equations numerically, see e.g. 
[ISl [2TI [3] . So far, the numerical evidence for a finite time blow-up is not yet conclusive 

pinsum. 

Global regularity results for the 3D Navier-Stokes equations have been obtained 
using energy estimates under some smallness assumption on the initial data [26 | l29 | 125] . 
It is well known that the convection term does not play an essential role in energy 
estimates due to the incompressibility of the velocity field. Thus, it is not clear how 
convection may contribute to global well-posedness of the 3D incompressible Navier- 
Stokes equations. 

In [21] , Hou and Li investigated the stabilizing effect of convection via an exact ID 
model. They found that the convection term plays an essential role in canceling the 
destablizing vortex stretching terms in this ID model. This observation enabled them 
to obtain a pointwise estimate via a Liapunov function which controls the dynamic 
growth of the derivative of vorticity. Motivated by the work of [2T], Hou and Lei 
further investigated the role of convection by constructing the following 3D model of 
the axisymmetric Navier-Stokes equations with swirl [22J: 

8 t u = u(d 2 + \d r + d 2 )u + 2ud z i>, 

d t u = u(d 2 + \d r + d 2 )u + d z (u 2 ), (1.1) 

-{d 2 r +ld r + dl)i, = uj. 

When v = 0, we refer to the above model as the 3D inviscid model. This model 
derived by using a reformulated Navier-Stokes equations in terms of a set of new 
variables (u,u,i/j) = (u e ,u e ,^p e )/r, where r = \Jx\ + x 2 ,, u e is the angular velocity, 
ui e the angular vorticity, and ip e the angular stream function, see [2TT [22] . The only 
difference between this 3D model and the reformulated Navier-Stokes equations is that 
we neglect the convection term in the model. This new 3D model shares several well 
known properties of the full 3D Euler or Navier-Stokes equations [22]. These include 
an energy identity for smooth solutions, an artificial incompressible constraint, a non- 
blowup criterion of Beale-Kato-Majda type [2], a non-blowup criterion of Prodi-Serrin 
type [301 E] , and a partial regularity result for the model [23J , which is an analogue of 
the Caffarelli-Kohn-Nirenberg theory [U [28] for the full Navier-Stokes equations. 

Despite the striking similarity at the theoretical level, this 3D model seems to have 
a very different behavior from that of the Euler or Navier-Stokes equations. Numerical 
study in [22] seems to suggest that the model develop a potential finite time singularity 
from smooth initial data with finite energy. However, the mechanism that leads to the 
singular behavior of the 3D model seems to be destroyed when the convection term is 
added back to the model. In a recent paper [23], Hou, Shi and Wang proved rigorously 
that the inviscid model can indeed develop a finite time singularity for a class of smooth 
initial data with some Neumann-Robin boundary condition. Moreover, they proved a 
global well-posedness result for a class of small smooth initial data. The analysis in [23] 



was carried out for a rectangular domain, which does not apply to the axisymmetric 
domain considered in this paper. 

In this paper, we prove that the 3D inviscid model with some appropriate boundary 
conditions of Neumann-Robin type can develop a finite time singularity starting from 
smooth initial data on a bounded or an unbounded exterior axisymmetric domain. 
Moreover, we obtain a global well-posedness result for a class of large smooth initial 
data under some Dirichlet boundary condition. 

Our analysis in this paper is similar in spirit to that of [23]. However, there are 
several new ingredients in this work. The first one is that our blowup result applies 
to a bounded or unbounded exterior domain in the axisymmetric geometry. The local 
well-posedness of the 3D model with a boundary condition of the Neumann-Robin 
type is more complicated than the case considered in [24J and involves the use of the 
modified Bessel function. The second ingredient is to construct a special positive test 
function that satisfies several requirements. In the case of a bounded domain, we use 
the first eigenfunction of the Laplacian operator with homogeneous Dirichlet boundary 
condition. In the case of an unbounded exterior domain, a positive and decaying 
eigenfunction of the Laplacian operator does not exist. This requires us to relax the 
constraints on the test function. The third ingredient is an improved estimate which 
enables us to establish the global regularity of the 3D inviscid model for a class of 
initial boundary value problem with large initial data. 

The paper is organized as follows: In Section 2 we will state our main results and 
present our main ideas of their proofs. Section 3 is devoted to proving the finite time 
singularity of the 3D inviscid model. In Section 4 we present our global well-posedness 
result. Finally, we prove the local well-posedness of the initial boundary value problem 
in Section 5. 

2 Main Results 

In this section we set up our problem and state our main results. First of all, let 
us recall that the 3D model ( ll.lj) is formulated in terms of a set of new variables 
u = — , u; = — , tp — — . Due to this change of variables, the original three-dimensional 
Laplacian operator for if) 6 has been changed into a five-dimensional Laplacian operator 
for ip in the axisymmetric cylindrical coordinate in ( 11. ip . Thus we can reformulate our 
3D inviscid model in the three-dimensional axisymmetric cylindrical domain 



£>(7i, 7 2 ) = {{x x ,x 2 , z) G M 3 : 71 < r = \jx\ + x\ < 72; < z < 1} 

as a five- dimensional problem in the axisymmetric cylindrical domain 

4 
^(7i,72) = {(si,-" ,x±,z)eM 5 :-/ 1 <r= {^xrf <7 2 ;0<z<l}. 



It is much more convenient to perform the well-posedness and the finite time blow-up 
analysis for our 3D model in the five-dimensional setting. In the remaining part of the 
paper, we will carry out our analysis in this axisymmetric five- dimensional domain. 
Denote 

f k r e - kcosh ^ cosh 2 (6)d6 A 

Now we are ready to state the first result of this paper which is concerned with the local 
well-posedness of classical solutions to the 3D inviscid model with a Neumann-Robin 
type boundary condition on the exterior domain Q(l, oo). 

Theorem 2.1. Let S exter i or be defined in ( 12.1 j) and (3 G ^exterior- Consider the initial- 
boundary value problem of the 3D model (II. ip with the initial data 

u G # 3 (ft(l, oo)), ip G H\n(l, oo)), (2.2) 

Uq > for z 7^ and z^l, Wo|«=o = u o\z=i = (2.3) 

and the Neumann-Robin type boundary condition 

(VV + W)|r=l = 0, V*l*=0 = *Pz\z=l = 0. (2.4) 

Then there exists a unique smooth solution (u, ip) to the 3D inviscid model with the 
initial data ( I2.2p -( l2~3j) and the boundary condition ( 12.41) on [0,T) for some T > 0. 
Moreover, we have 



u 



eC([0,T),# 3 (fi(l,oc))), ^GC([0,T),iJ 4 (fi(l,oo))). (2.5) 



The proof of Theorem 12.11 relies on an important property of the elliptic operator 
with the mixed Neumann-Robin type boundary condition. Note that the boundary 
condition in (12.41) is not the third type. To recover the standard elliptic regularity 
estimate, we need to study the spectral property of the differential operator which 
requires that (3 G ^exterior- We remark that the local well-posedness analysis for (II .ip 
with a mixed Dirichlet-Robin boundary condition has been established in [23] where 
the non-standard boundary condition is imposed on the axial direction. Here our 
non-standard boundary condition is imposed on the radial direction. The local well- 
posedness in such case involves the use of modified Bessel function and is more involved. 
The proof of Theorem 12. II will be deferred to Section 5. 

Next, we state the finite time blowup result of the exterior problem (II. ip . (l2.2IUI2.3p 
and (El). 

Theorem 2.2. Suppose that all the assumptions in Theorem \2.1\ are satisfied. Let 
f3 > 2 + 2^1 + ^, a = f and 



(r, z) = 9{r) sinM, 9{r) = e~ ar \ $(r, z) = cj> rr + -<j> r + <f> gg . (2.6) 
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Then the solution of the 3D inviscid model will develop a finite time singularity in the 
H 3 norm provided that 



M /*00 /»! /*00 

/ / $ \og(u 2 )r 3 drdz > 0, / / $ip 0z r 3 drdz > (2.7) 

Jo J\ Jo J\ 



and 

r 1 r°° \2 iq / r 1 r°° x 3 

/ / d z ^ ®r 3 drdz) > — ( / / (log u 2 )<5>r 3 drdz) , (2.8) 

where Cq is a positive constant defined in (I3.14p . 

The proof of Theorem 12.21 is in spirit similar to that of Theorem 1.1 in [24J. The 
main ingredient here is to construct a special positive test function which meets our 
requirements on the unbounded domain case. In [21], the authors can take the product 
of sine functions in the x — y domain as a testing function since the domain is a 
rectangular domain. In the case of a unbounded exterior domain in the axisymmetric 
geometry, a positive and decaying eigenfunction does not seem to exist. This requires 
us to relax the constraints on the test function. We will present the details of the proof 
of Theorem 12.21 in Section 3. 

Remark 2.3. We remark that similar results in Theorem \2.1\ and \2. 6 A are also true 
in the bounded domain case Q(0, 1). We will not present this result here. A more 
interesting result would be to obtain a similar result for the case of Q(0, 00) with an 
initial boundary condition whose energy is conserved in time. This would require a 
different technique. We will report it in a forthcoming paper. 

In Theorem 1.1 of [21], the authors proved the finite time singularity of the 3D 
inviscid model in a rectangular domain. Their Robin boundary condition is imposed 
on the axial boundary z = 1. Our next theorem extends their result to the case of a 
bounded axisymmetric domain Q(0, 1). 

Let (Afc, 0k{r)) be the eigenvalue-eigenfunction pair of the following Dirichlet eigen- 
value problem: 

X k 9 k , re [0,1), 

Define 

interior = {P > | f3 ± Afc, f3 ^ ^^ + e ^ for all k = 1, 2, ■ • • }, (2.10) 




gV A fe — g~ V A fc 



and 

0(r, z) = — d x {r). (2.11) 



e -a(z-l) _|_ £ a(z-l) 



Moreover, we assume that a, (3 satisfy 

We can prove the following local-well-posedness and finite time blow-up results. 

Theorem 2.4. (A) Local well-posedness. Let ^interior be defined in (j2.10p and (3 G 
interior- Consider the initial- boundary value problem of the 3D model (II. ip w£/i t/ie 
initial data 

u eH 3 (tt(Q,l)), ^ Gif 4 (fi(0,l)), «oU=o = «o|*=i = (2.13) 

and the Dirichlet-Robin type boundary condition 

iP\ r=1 = i,\ z=l = 0, ip z + M\ z=0 = 0. (2.14) 

TTien t/iere exists a unique smooth solution (u, ip) to the 3D inviscid model with the ini- 
tial data (I2.13P and the boundary condition (I2.14p on [0, T) for some T > 0. Moreover, 
we have 

ueC([0,T),H 3 (n(0,l))), ^GC([0,T),tf 4 (fi(0,l))). (2.15) 

(B) Finite time blow-up. Suppose (I2.12p is satisfied and u\ > for < z < 1. //% 
and ^o satisfy 



(logw o )0r drdz > 0, / / ijJo z 4>r drdz > 0, 
(I jo Jo Jo 



and 



( / / d z tpo4>r 3 drdz) > —( / (logWo)0r 3 drd2: 

for some absolute constant c\ > 0, then the solution of the 3D inviscid model will 
develop a finite time singularity in the H 3 norm. 

The local well-posedness of the initial-boundary value problem for 3D model in 
Theorem 12.41 can be carried out using almost the same argument as that in [23J by 
replacing the first eigenfunction sin nxi sin 11x2 sin irx^ sin nx^ of the rectangular domain 
by the first eigenfunction 4>i(x) of the unit ball in M 4 . The finite time blow-up result 
can be proved in exactly the same way as in [23]. Since the analysis is essentially the 
same as that in [24], we will omit the proof of Theorem 12.41 in this paper. 

The last theorem extends the global existence result in Theorem 6.1 in [23] to 
the axisymmetric domain Q(0, 1). The most interesting aspect of this result is that 
the initial condition can be made as large as we wish in the Sobolev space. This is 
achieved by effectively imposing a sufficiently large negative boundary condition for 
dzi>\dri(o,i) — —M with M being a large positive constant. 



Theorem 2.5. Let M > an arbitrarily large constant and s > 3. Assume that 
u G H s (n(0, 1)) ; tp G H s+1 (n(0, 1)) with uo(r, 0) = w (r, 1) = and Vo|r=i = -M«, 
<9 2 "00 |z=o = d.zV'ol^i = —M. Then the 3D inviscid model with the initial and boundary 
data 

«(O,.)=«o(0, V(o,-)=M0> (2-16) 

and 

4>{t, 1, 0) = -Mz, <9^(t, r, 0) = d z 4>{t, r, 1) = -M (2.17) 

admits a unique global smooth solution (u,ip) with u G C([0, oo); if s (fi(0, 1))), ^ G 
C([0, c»); # s+1 (fi(0, 1))) provided that 

||V9^o|U.-i < M/(8C 2 S ), \\u 2 \\ H s < M 2 /(AC 3 S ), (2.18) 

where C s is an absolute positive constant depending on s and Q only. Moreover, we 

have 

\\„,u\ 2 \\ <r n IU, 2 II „-2Mt icy 1Q \ 

and 

C 
1^9^11^-1(0(0,1)) < HV^VollH'-ifnco.i)) + ^ttII m oII^ s (^(o,i))- (2.20) 

We will present the proof of Theorem 12.51 in Section 4. 

3 Finite Time Singularity in the Exterior Domain 



In this section, we present the proof of Theorem 12. 2\ which shows that the 3D inviscid 
model develops a finite time singularity in the exterior domain. 

Proof of Theorem 12.21 We will prove Theorem 12.21 by contradiction, as in [21]. By 
the local well-posedness result, we know that the initial boundary value problem of the 
3D inviscid model has a unique solution u G H 3 (Q(1, oo)) and ip G H A (Q(l,oo)) for 
< t < T for some T > 0. Let T& be the largest time for which our 3D inviscid model 
has a unique smooth solution u G H 3 (Q) and ip G H 4 (Q) for < t < T&. We will prove 
that Tfe < oo. Suppose that T& = oo. Then we have 



« 



eC([0,oo),H 3 {Q{l,oo))), ^eC([0,oo),# 4 (ft(l,oo))). (3.1) 



We will prove that (13. ip would lead to a contradiction. 

Define and $ as in Theorem 12.21 By a straightforward calculation, we have 

$(r, z) = [Aa 2 r 2 - (8a + vr 2 )]0(r, z). (3.2) 



n-2 



For a > 1 + \/l + ^-, it is easy to get that 

4aV - (8a + n 2 ) > Aa 2 - (8a + ix 2 ) > for r > 1. 



Consequently, we have 

$(r,z)>0 for 1 <r,0 < z < 1. (3.3) 

On the other hand, due to the boundary condition on the initial data Uq in (I2.3p . one 
can use the first equation in (II -ip with v = to solve u: 

u 2 (t,r, z) = ul(r, z) exp{4 / d z ip(s, r, z)ds} for z ^ 0, 1. 

Jo 

Hence, by continuity, as long as the solution is smooth, one also has 

u(t, r, 0) = u(t, r, 1) = 0. (3.4) 

Multiplying the second equation in (11.11) with v = by <p z and integrating over 
fi(l, oo), we have 

1 /*00 q /*1 /*00 

/ (d 2 + -<9 r + d 2 z )ip t (f) z r 3 drdz = / / d z u 2 (j) z r 3 drdz. (3.5) 

'0 A r A A 

On one hand, using H3.4I) and performing integration by parts, we get 

pi poo pi poo /»1 /*oo 

/ / d z u 2 (fi z r 3 drdz — — / u 2 d 2 z (fir 3 drdz = n 2 / u 2 (f>r 3 drdz. (3.6) 

Jo A Jo Ji Jo Ji 

On the other hand, using the boundary condition (12.41) and performing integration by 
parts, we obtain by noting that 8 r (l) + (36(1) = that 

pi poo q 

- / / (d 2 r + -d r + d 2 )^ z r 3 drdz (3.7) 

d f 1 d f l f°° 3 . 

= ~n {^r^z - i)d 2 rz 4>)\ r= idz - — / ijj(d 2 + -d r )<p z r 3 drdz 
dt J dt J J 1 r 

d f l Z" 00 
+ — / / ip z d 2 (j)r 3 drdz 
dt Jo J i 

d f 1 d f 1 f°° 

= n 9(1)— (cosTTz)(ip r + (3ip)\ r= idz + — / / ip z $r 3 drdz 
dt Jo dt Jo J x 

d f 1 r°° 

= — / / i/j z ^>r 3 drdz. (3.8) 

dt Jo Ji 

Combining (13. 7p with ( 13 . 5 [) - ( T3T6]) gives 

j pi poo pi poo 

— ip z $r 3 drdz = tt 2 / / u 2 <j)r 3 drdz. (3.9) 

at Jo Ji Jo Ji 



Multiplying the first equation in ( II. ip with v = by $ and integrating on Q(l, oo) 
yield 

j pi poo pi poo 

— / / (logu 2 )$r 3 drdz = A / / ip z $r 3 drdz. (3.10) 

Since JI Jj 00 d z ip $>r 3 drdz > 0, we conclude from (13 .9p that 

/ ip z $r 3 drdz>0 for all t > 0. (3.11) 

It follows form (I3.10p - (l3.1ip and the condition J^ J^° (log u^)^r 3 drdz > that 

/ / (logu 2 )$r 3 drdz > for all £ > 0. (3.12) 

Consequently, we obtain by using (I3.2p and (" 13. 12[) that 

( / / (logu 2 )$r 3 drdz\ < ( (log + u 2 )$r 3 drdz\ (3.13) 

< 4( / / \u\$r 3 drdz) 

pi pOO pi pOO jTj 

<4 / / \u\ 2 4>r 3 drdz / / (—) 2 <j)r 3 drdz 
Jo Ji Jo Ji <P 

8c 7r 2 f 1 f°° . l2j „, , 
< / / \u\ wr drdz, 

~ 3 Jo Ji l 
where c is defined by 

■1 POD 



3 
Co 



2tt 2 

Using (EU), ( 13TT0]) and (13TT3]) . we get 
d 2 ^ '°° 



/ / (4aV - (8a + ir 2 )) 2 (j)r 3 drdz < oo. (3.14) 



dt ./ o j i 



Let 



Then (13. 15[) implies 



/•I /"OO 

2 )$r 3 drdz = An 2 / / u 2 (f)r 3 drdz 


(3.15) 


3 / f 1 f°° . \2 

> ( / / (log u 2 )$r 3 drdz) . 

2c V Jo J\ ' 




rl /*oo 

'(*)=/ / (log u 2 )$r 3 drck. 


(3.16) 


y"(t) > i-y(t). 


(3.17) 



It follows from (13~T0>(l3TTTj) that Y'(t) > 0. Multiplying (IXTTj) by Y'(t) and integrating 
in time from to t, we get 

(Y'(t)Y >(Y'(0)) 2 --Y(0f+ -Y(t) 3 
\ J \ J c c 

= Uf f d z ip ^r 3 drd z y -—( J J (loguD^drdzY + —Y(t) 3 . 

The condition ( I2.8P implies 

(Y'(t)Y > -Y(t) 3 . (3.18) 

V / c 

Since Y'(t) > 0, it is easy to solve (13. 18[) to conclude that 

/ / (\ogu 2 )$r 3 drdz> Jo Jl V 6 o; - (3.19) 

Jo Jl (2 v ^-t v // 1 / 1 oo (log^)$ r 3^ 

Note that (Q gives 

$(r, z) = [Aa 2 r 2 - (8a + 7r 2 )]0(r, z) = [Aa 2 r 2 - (8a + n 2 )]e- ar2 sin(Trz). (3.20) 

Since r > 1, it is easy to show that $ satisfies 

< $(r, 2) < 4a 2 r 2 e" ar2 < 4aV°, r > 1. (3.21) 

Thus, we have 

pi poo pi poo 

j j (\ogu 2 )$r 3 drdz < (\og + u 2 )$r 3 drdz (3.22) 

Jo Ji Jo Ji 

/l poo 
/ u 2 r 3 drdz. 

Combining (13. 19[) with (I3.22J) gives 

9 „ f 1 f°° 9 ^ , , 4c o (r! fi°° (log ul)$r 3 drdz 

Aa 2 e~ a / / u 2 r 3 drdz > Jo Jl v 6 o; _ . (3.23) 

Jo Jl (2^-t^ l J™{\ogu 2 )$rHrdz) 

On one hand, we have j f ± u 2 r 3 drdz < oo for all times by (j3.ip . On the other 
hand, the right hand side of ( I3.23I) will blow up as t — > , -, = . This is 

clearly a contradiction. This contradiction implies that the assumption ( 13. ip can not 
be valid and the solution must blow up in a finite time in the H 3 norm no later than 
T* = . -. v co =. This completes the proof of Theorem [2] 

10 



4 Global Regularity with Large Data 

This section is devoted to proving Theorem 12.51 

Proof of Theorem 12.51 First of all, the local well-posedness of the initial boundary 
value problem can be established by using an argument similar to that of [24J . Now let 
us assume that (u, ip) is a local smooth solution satisfying boundary condition (" 12 . 1 T[) 
such that u G C([0, T); H s (tt(0, 1))), ip G C([0, T); H s+1 (tt(0, 1))) for some T > 0. We 
are going to show that T = oo. 
Denote by 

v = —d z ip, u = u 2 . 

It is easy to see that 

v G C([0, T); iT(fi(0, 1))), « G C([0, T); tf* (fi(0, 1))) • 
Differentiating the second equation in ( 11. ip with u = with respect to 2, we get 

(9 r 2 + -9 r + 9 Z 2 )^ = 9 2 2 m. 
r 

Note that t?t^T|©rj(o,i) — 0. We define g = A^" 1 / as the solution of the Laplacian equation 
with the homogeneous Dirichlet boundary condition: 

4 
C£dl+d 2 z )g = f, 010,(0,1) = 0, V/a 2 (U(0,l)). 
t=i 

Then, we can reformulate the 3D inviscid model as follows: 

u t = —Auv, 
v t = A^d 2 u. 

Note that the boundary condition in (12. 17ft implies that 

v\aci(p,i) = M - 

If we further denote 

V = v-M, (4.1) 

we obtain an equivalent system for u and v. 

u t = -AMu - Auv, (4.2) 

iH = A-'d^u, (4.3) 

v\an(o,i) = 0. (4.4) 

Recall the following well-known Sobolev inequality [15]. Let u, v G H S (Q) with 
s > n/2 (n is the dimension of Q). We have 

\\uv\\ H °(n) < C s \\u\\ H °(n)\\v\\H°(n)- (4.5) 

11 



We will also use the Poincare inequality 

IMIif(n) < C s \\Vv\\ H s-i (n) , (4.6) 

if v \dCi(o,i) — 0. Here C s is an absolute positive constant depending on s and Q only. 
Taking the H s norm to the both sides of ( 14. 2 p and using (J4.5p ~ (j4.6p . we obtain 

-tJu\\h° < -4M\\u\\ H . + 4C a \\u\\ H .\\v\\ H . (4.7) 

at 

< -4M||u||ff. +4:Cs\\u\\ H s\\Vv\\ H s-i 

< (-4M + 4C?||Vv|| JT .-i)||u|| H ., 

where we have used (14.41) . 

Next, we apply V to the both sides of ( 14. 3ft and take the i? -norm. We get 

j^VvWhs-! < ||VA^u|| H .-i (4.8) 

< C s \\d 2 z u\\ H s-2<C s \\u\\ H s, 

where we have used the standard elliptic regularity estimate for Aj 1 . 

By the local well-posedness result and the assumption on the initial condition ( 12. 18ft . 
we know that there exists a positive T > such that we have 

M 
IIV^Hh,-^— , 0<t<T. (4.9) 

Let T* be the largest time such that ( 14. 9 p holds. We will prove that T* = oo. Suppose 
T* < oo. Substituting (gU) into M~T} . we obtain 

— \\u\\ H s < -2M\\u\\ H s, 0<t<T*, (4.10) 

which implies 

INU* < \\u \\h^- 2M \ 0<t<T*. (4.11) 

Now substituting (T4TTTT) into P~%|) yields 

||Vu||h.-i < \\Vv \\ H s-i+C s \\u \\ H s e~ 2Ms ds (4.12) 

Jo 

C 

< UVuollff-i + 2j^||wo||h. 

M M 

where we have used the condition ( I2.18p . Since ||Vt>||# s -i < -^ < -^ for t < T* 

and ||w||#« < ||iTo||// s e _2MT * < ||wo||.ffa for t < T*, we can use our a priori estimates 
( |4\TP -( I4~8]) to show that (TC9]) remains valid for < t < T* + 5 for some 5 > 0. This 
contradicts with the assumption that T* is the largest time such that ( 14. 9 p holds. This 
contradiction implies that T* = oo and estimates ( l4.lip -f T4.12p remain valid for all 
times. This completes the proof of Theorem 12.51 
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5 Local Well-Posedness in the Exterior Domain 

The local well-posedness theory of the initial-boundary value problem of the 3D model 
(II. ip with ( 12. 2 p and (12. 4p is based on the following lemma 15.11 Once the lemma is 
proved, the proof can be carried out in exactly the same way as the local well-posedness 
analysis presented in [24] . Hence, we will only present the proof of the lemma [57T1 in a 
slightly modified domain with 1 < r < oo and < z < n. 

Lemma 5.1. For any given u G H S ~ 2 (Q(1, oo)) with s > 2, there exists a unique 
solution ip G H S (Q(1, oo)) to the boundary value problems 

(<9 2 + fd r + dfji) = u, 1< r, < z < vr, 

tfr + fti>\r=l = 0, 1p z \ z =Q = ^z\z=tt = 0, 

where (3 G Exterior is a constant. Furthermore, we have the following estimate: 

\\iP\\h»(ciq.,oo)) < C|MlH s - 2 (n(i,oo)), (5.2) 

where C is an absolute positive constant. 

Proof. Let us decompose ip = ip^ l > + ip^ 2 \ where ip^ is the solution of the elliptic 
equation with the following mixed Dirichlet-Neumann boundary condition: 



-(<9 2 + f«9 r + <9 2 )^ (1) =w, infi(l,oo), 
i>P |*=o = i/>P\s=* = 0, ^ (1) |r=i = 0, ^ (1) -> as r ->■ oo. 



(5.3) 



Then ?/>( 2 ) satisfies 



(«9 2 + f,9 r + ,9 2 )^ 2 ) = 0, infi(l,oo), 

[4 2 ) + ^(2)]| r=1 = _^)| r=1; (5.4) 

,4 2) | z= o = 4 2) | 2 =. = 0. 

The standard elliptic theory gives that ip^ 1 ' G H s (Q(l, oo)) and H^ ||if s (n(i,oo)) < 
IM|ff s - 2 (n(i,oo))- It remains to show that the Neumann-Robin problem (15.41) is well- 
posed. 

Let us perform the cosine transform to ( 15. 4ft with respect to z variable, which gives 

'-(<9 2 + f<9 r )^ + A; ; V( 2 ) = 0, 



where the cosine transform ^ is defined as: 

2 



-0( 2 )(r, k) = — / ip^ '(r, z) cos(kz)dz. (5-6) 

7T 
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Let ?// 2 )(r, k) = -f(rk, k). Then (15.51) can be written in terms of / as 

k2 ^ k ) + k -ifW ~ (k 2 + ±)f{rk) = 0, 

which is equivalent to the following modified Bessel equation 

r 2 f"(r, k) + rf'(r, k) - (1 + r 2 )f(r, k) = 0. 

The general solution of the modified Bessel equation with the decaying property as 
r — y oo is given by modified Bessel function: 

f(r,k) = C(k)K(r), 

where K(r) is the modified Bessel function and can be represented in an integration 
form as follows 

/•oo 

K{r) = / e - rcosh(e) cosh(0)d0, K(r) -y as r ->■ oo. 
Jo 

In fact, we have < K{r) < C2r~ 1 ^ 2 e~ r as r — y oo for some positive constant ci- Thus, 
the general solution of ip^(k,r) is given by 

^)(r,k) = C ^-K{rk). (5.7) 

r 

The boundary condition of ip^ 2 ' in ( 15. 5ft can be used to determine C(fc) as follows: 

Note that C(k) in ( 15.81) is well-defined for (3 G Exterior- The remaining part of the proof 
is exactly the same as in [23]. We can show that ip^ will have the same regularity 
property as ip^\ This proves our lemma. We will omit the details here. □ 
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